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Abstract 

This paper presents a new minimum-weight 
trellis-based soft-decision iterative decoding al- 
gorithm for binary linear block codes. Simu- 
lation results for the RM(64,22), EBCH(64,24), 
RM(64,42) and EBCH(64,45) codes show that the 
proposed decoding algorithm achieves practically 
(or near) optimal error performance with signif- 
icant reduction in decoding computational com- 
plexity. The average number of search iterations 
is also small even for low signal- to-noise ratio. 

1 Introduction 

Recently Moorthy et. al. [2] have proposed a 
zero-and-minimum-weight subtrellis-based itera- 
tive decoding scheme for binary linear block codes 
to achieve a very good trade-off between error 
performance and decoding complexity. In the 
scheme, all the candidate codewords are gener- 
ated by an algebraic decoder based on a set of 
test error patterns proposed by Chase[l]. The 
zero-and-minimum-weight trellis search around 
the current best candidate codeword c is per- 
formed at most once, only if (i) a sufficient condi- 
tion that the optimal solution is within the min- 
imum distance from c holds or (ii) all the test 
error patterns have been exhausted and no can- 
didate codeword satisfies the sufficient condition 
for optimality. 

For the proposed decoding algorithm in this 
paper, preliminarily presented in [6], the initial 
candidate codeword is generated by a simple de- 
coder, the zero-th order or the first order decod- 
ing proposed in [4]. The subsequent candidate 
codewords (if needed) are generated by a chain of 
minimum- weight trellis searches. This minimum- 
weight trellis around a candidate codeword c con- 
sists of only the codewords in code C that are at 
the minimum distance from c, but does not in- 
clude c. The decoding iteration stops whenever a 
candidate codeword is found to satisfy a sufficient 
condition for optimality or the latest minimum- 
weight trellis search results in a repetition of a 
previously generated candidate codeword. Let 
this decoding algorithm be denoted Algorithm I- 


w\. The decoding process terminates faster than 
the Moorthy et. ah algorithm. Furthermore, the 
use of minimum- weight trellis search considerably 
reduces the possibility of being trapped into a lo- 
cal optimum. As a result, it achieves better 
error performance than the Moothy et. al. algo- 
rithm. 

A necessary condition for Algorithm I-uq to 
achieve good error performance is that the min- 
imum weight codewords span the entire code. 
Reed-Muller(RM) codes satisfy this condition. 
Simulation results for the RM(64,22)(the (64,22) 
RM code), RM(64,42) and the EBCH(64,45)(the 
extended (64,45) BCH code) codes show that the 
proposed decoding algorithm practically achieves 
optimum MLD performance even in the range of 
relatively low SNR. The EBCH(64,24) code is an 
example for which the above necessary condition 
does not hold. For this code, the first and second 
minimum weight codewords span the entire code. 
In this case iterative decoding algorithm based 
on the first and second minimum-weight trellis 
search, denoted Algorithm \-w\-W 2 , is used. 

We also propose another approach to overcome 
the problem. Let Co be a linear subcode of C and 
assume that the minimum weight codewords of Co 
span Co. The decoding scheme is a combination 
of: (1) the iterative search using the minimum- 
weight trellis for Co around the latest candidate 
codeword, and (2) a procedure for moving from 
the coset of Co in C, containing the current candi- 
date codeword, to another coset which is likely to 
contain the optimal solution. Simulation results 
for the EBCH(64,24) code show that this scheme 
achieves better error performance than Algorithm 

2 Sufficient Conditions for Op- 
timality 

Suppose a binary (TV, K) linear block code C is 
used for error control over the AWGN channel 
using BPSK signaling. Let z = (zi,Z 2 » • • • > z n) be 
the binary hard-decision sequence obtained from 
the received sequence r = (r i , 7*2 , . . . , r jv ) . 

Let V N denote the vector space of all binary 
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//-tuples. For an //-tuple u = (ui,U 2 - ■ • u^) E 
V N , let L(u) be defined as follows: 

L(u) = ^2 M- (!) 

{i-.ui^zi and i<i<N} 

L(u) is called the correlation discrepancy of u 
with respect to z, and the smaller L(u) is, the 
larger the correlation between u and r is. For 
u and v E V^, u is said to be better than v 
if L(u) < L(v). For a nonempty subset X of Vn 
and a positive integer h , let h! denote min{|X|,/i} 
and let best^X denote the set of the h! best n- 
tuples in X, that is, for any u E best^X and 
v E X - bestfcX, L(u) < L(v). The best n-tuple 
in X will be denoted bestX. 

Let (1h{u, v) denote the Hamming distance be- 
tween two N-tuples, u and v. For Ui, u 2 , 

. Uk E Vn and positive integers di, d 2 > • • 
d/t, let VN{ui,di\U 2 ,d 2 \ . . .\Uh,dh) be defined 
as the set: {tt E Vn : d B (u y Ui) > d x for 

1 < i < h}, and let L(ui,di; ti 2 *d 2 ; . . . ; 
be defined as the minimum of L(u) over tx E 
Vjv( < wi,di;u 2? d 2 ;...;u^ j d/ l )* 

Then we have the following early termination 
condition of an iterative decoding algorithm with- 
out degrading the error performance. 

Lemma: At a stage of an iterative decoding al- 

gorithm for a block code B(C itself or a coset of 
a linear subcode of C), let GC denote the set of 
those candidate codewords which have been gen- 
erated already, and let Ubest denote the best of 
GC. Suppose that for U\, u 2 > • • •> u h £ GC and 
positive integers di, d 2 , . d />, u be st is the best 
of ujL 0 {u E B : d H {u,Ui) < dj. If ti bes t satisfies 
the following condition, 

L(ubest) < L(ui, di; ti 2 ,d 2 ; . . . ; (2) 

then Ubest is the optimal MLD solution in B. AA 
Expressions for evaluating the right-hand side 
of (2) for h = 1,2 and 3 have been derived [3]. 
These bounds are used in our proposed iterative 
decoding algorithm for early termination condi- 
tions. 

3 Decoding Procedure 

Let C be a binary linear (N y K) code with weight 
profile W = {0,uq,tu 2 ,.. .}, and let Co be a bi- 
nary linear (X, Ko) subcode of C with weight pro- 
file Wo = {0,tuoi>wo2> . • *}> where w\ is the mini- 
mum weight of C, u>oi is that of Co and w x < WQi, 
As a special case, Co may be C itself. 

3.1 Minimum- weight Subtrellis Search 

for a Coset 

For B E C/Co and u 6 B, minimum- weight 
search (u, B) denotes a search procedure for find- 


ing a codeword in B , denoted which has 

the least correlation discrepancy with respect to z 
among all codewords in B at the minimum Ham- 
ming distance tuoi from u. That is, 

<Pb{u) = best{v E B : d B {v y u) ~ t^oi}* (3) 

If Co is spanned by the set of minimum weight 
codewords of C, that is, Co = C or C is an 
UEP(unequal error protection) code, then woi = 
W\ and 

<Pb(u) = <Pc{u) — best{u E C : dn{v, u) = w,}. 

( 4 ) 

This search procedure is implemented by using 
the minimum- weight subtrellis of Co around u. 
This minimum-weight subtrellis is sparsely con- 
nected and much simpler than the full trellis of 
the code[2, 6]. 

Iterative minimum-weight search (u, B) is to 
generate a sequence of candidate codewords, 

<p B (u), ( Pb{ ( Pb{ u ))i •**> until a certain termina- 
tion condition holds. It is shown in [6] that 

L(<^+ 2 (u)) < L(<p' B (u)) y for i > 0, (5) 

where <p%{u) = u and y X Q l {u) = ^(^b( u )) f° r 
i > 0, and that if L(<p l g 2 (u)) < L(<p l B (ii)) for 
0 < i < J, then <p l B {u) with 1 < i < I are all 
different. If j is the smallest index such that 

L« 2 (u)) = L(<p 3 b (u)), (6) 

then 

min{L(<^(u)) : 0 < i < j) = 

min {L(<p 3 q 2 {u)),L{^ 1 {u))}. (7) 

The condition (6), denoted Cond#, is used as one 
of termination conditions for Iterative minimum- 
weight search {u,B) to avoid repetition. 

3.2 Decoding Procedure for C 

Suppose the set of codewords of weight tuoi in Co 
spans Co and K - Ko is not large. We propose 
a new decoding procedure for C which consists 
of iterated minimum- weight searches in a coset 
and coset shiftings. Two early termination condi- 
tions, Condc for the entire procedure and Conda 
in the subprocedure for a coset B E C/Co, are 
used besides the termination condition Cond# in 
the subprocedure for a coset. 

Condc is a sufficient condition [3] that the best 
candidate codeword, denoted tibest* i n the set of 
those candidate codewords which have been gen- 
erated already, denoted GC, is optimum based on 
best^GC and the weight profile W of C, where h 
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is a specified small integer. From the Lemma, 
Condc is defined as 

L(ubest) < L{ui,di-,U2,d2\..-]Ui,di), (8) 

where l = min{/i, |GC|} and for 1 < i < /, Ui G 
besthGC and if <pc(«i) € GC, then d t = w 2 and 
otherwise, d x — w\. 

Cond# is a sufficient condition that there re- 
main no codewords in B better than Ubest* 
This condition is also based on best^GC# where 
GC# = GC n B and the weight profile Wo of Co 
which is the same as the distance profile of any 
coset of C/Co- From the Lemma, Cond# is de- 
fined as 

^( w best) 5: Li u U ^li u 2> * * • 5 u Ib > 

where {# = min{/i, |GC#|} and for 1 < i < Ib, 
U{ £ best hGCb and if <pb{ u \) € GC#, then d* = 
tt;o2 and otherwise, d{ = u>oi- For instance, if the 
minimum or the second minimum weight of C 0 
is greater than that of C, then Cond# is more 
effective than Condc only. 

In the procedure, global variables GCh and 
GCb^h are used besides GC and ^best- GCh de- 
notes the current value of best^GC and GC#^ 
denotes the current value of best^GC#. For 
B £ C/G 0 , let f{r,B) denote the initial candi- 
date codeword in B for a given received sequence. 

3.2,1 Decoding Algorithm II 
We assume that z ^ G. 

(Dl) (i) Generate f(r y B) for all B £ C/Co, and 
number the 2 K ~ K ° cosets of C/Co in the increas- 
ing order of correlation discrepancy, L(f (r , B)) - 
(ii) Initialize GC «— {f(r y B) : B £ C/Co}, 
Ubest 4- bestGC and GC h f- best^GC. If Condc 
holds, then output “Ubest an d stop. Otherwise, ini- 
tialize GC#, GCs,h 4— {/( r >F0} f° r & ^ C/Co, 
and perform Search-in(the first coset). 

(D2) Search-in (S): Execute Iterative minimum- 
weight search(/(r, B), B) together with updat- 
ing the global variables each time a new candi- 
date codeword is generated until either Cond#, 
Condc or Cond# holds. If Condc holds, then 
output Ubest and stop. Suppose that Cond# or 
Cond# holds. If all cosets have been exhausted, 
then output Ubest and s ^°Pj aQ d otherwise, call 
Search-in(the coset next to B). 

(D3) If either Cond c or Cond# holds for every 
coset in C/Co, the output is optimum. 

For the special case where C = Co, this decod- 
ing algorithm becomes Algorithm I-uq. 

3,3 Choice of the Initial Canidate 
Codeword /(r,B) for B £ C/Co 

For a given received sequence r = *• ,rjv), 

let Mk be the location set of the most reliable 


basis of the column space of a generator matrix of 
C, and let A;v_ # 0 the location set of the least 
reliable basis for the column space of a parity- 
check matrix of Co- Then it follows from Theorem 
1 in [5] that 

\M K nA N - Ko \ = K-K 0 . ( 10 ) 

For x = (x 1 ,x 2 ,...,xjv) € V N and a coset B G 
C/Co, a codeword u = (ui,u 2 , ■ ■ ■ ,un) G B sat- 
isfying the following condition is uniquely deter- 
mined: 

m = Xi, for all i G Mk ~ An-Ko- (H) 

Let g(x,B) denote the above codeword u in B. 

Then g(z,B) can be chosen as the initial candi- 
date codeword f(r,B) for B = C/Co, where z is 
the hard-decision binaxy vector obtained from the 
received sequence r. This g(z,B) for B = C/Cq 
is a simple generalization of the zero-th order de- 
coding proposed by Fossorier and Lin[4] to a coset 
of Go in C. Similarly, a generalization of their first 
order decoding can be used. 


4 Examples 

Example 1: Let C = EBCH(64,24) and Go — 
RM(64,22). Then w i = 16, w 2 = 18, w 0 \ = 16 
and to 0 2 = 24. Decoding Algorithm II, where 
early termination condition Condg is not used, 
has been simulated for this code. The simulation 
results are shown in Figures 1 and 2. For compar- 
ison, the simulation results of Algorithm I-uq-u>2 
for this code are also shown. Figure 1 shows the 
bit error probabilities. We see that Algorithm 
II practically achieves optimal error performance. 
Figure 2 shows the average numbers of opera- 
tions (addition and comparison) of Algorithm I- 
wi-w 2 and Algorithm II. The numbers depend on 
the complexities of subtrellises used in the sim- 
ulation. The minimum weight subtrellis of the 
RM(64,22) code used in Algorithm II and the first 
and second weight subtrellis of the EBCH(64,24) 
code used in Algorithm l-wi-w 2 are obtained sim- 
ply by purging the 4-section full trellis diagrams 
of the codes. The construction of better subtrellis 
is under study. 

We see that the average number of operations 
of Algorithm II can be reduced by using Cond# . 
Example 2: Let C = C 0 = EBCH(64,45) with 
Wl = 8 and w 2 = 10. For this case, the simula- 
tion results of Algorithm 1-uq, where the initial 
candidate codeword is provided by the first order 
decoding in [5], are shown in Figures 3 and 4. 

Simulation results of the RM(64,22) and 
RM(64,42) codes are shown in [6]. 
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Figure 1: Bit error probabilities for the 

EBCH(64,24). 



Figure 2: Average numbers of operations for 

EBCH(64,24). 
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